Dynamic hysteretic features of the Ising-type thin films 
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In order to elucidate the important characteristics of hysteretic response, such as types of the 
frequency dispersion curves, the decay of hysteresis loop area, the mechanism of domain nucleation 
and/or growth in the dynamic process, etc., a foUowup examination after the original work has 
been presented. The additive complementary treatment is essentially based on the results of recent 
preprint arXiv:1302.2727 concerning the effect of oscillatory perturbation on the surface enhance- 
ment phenomenon in magnetic thin films. Throughout the analysis, the best appropriate parameter 
values have been chosen since they would allow us to achieve striking results. The topology of the 
response for different layer indices in two regimes of modified surface exchange has been particularly 
emphasized. 
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I. INTRODUCTION 

Hysteresis is the signature of how the cooperative 
many-body system parameters respond dynamically to 
an alternating external magnetic field sweep. For the ki- 
netic Ising-type systems, the most familiar one is orig- 
inates as the ordinary magnetization (time-dependent 
magnetization) versus external field in the form of Lis- 
sajous curve. This ubiquitous phenomenon has been sub- 
ject of intensive interest due to a broad range of appli- 
cations. Hysteretic behavior is, however, a complicated 
process of a dynamic and nonlinear nature that elude 
serious treatment, both experimentally and theoretically 
since the area enclosed by the hysteresis loop (HL) is di- 
rectly proportional to the energy loss in a magnetization- 
demagnetization cycle. The remanent magnetization (in 
other words residual magnetization which is the magne- 
tization left behind in the system after an external mag- 
netic field is removed) and coercive field which means 
that the intensity of the external magnetic field needed 
to change the sign of the magnetization has been calcu- 
lated at several times in literature by benefiting from the 
hysteresis in order to understand and clarify the behavior 
of the system under consideration. HLs areal scaling law 
-which enables universality classification- has been mani- 
fested first by Steinmetz [1] empirically on his pioneering 
work dating back to the end of the nineteenth-century. 

Many efforts have been devoted to prediction and ex- 
perimental verification of the scaling behavior of the hys- 
teresis loop area (HLA) for thin films (for a brief review 
of HLA scaling results see Ref . [2] ) • The recent series of 
studies about characteristics of the hysteretic response 
in thin films (the topic also constitutes the essence of 
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our examination) have been also propounded by differ- 
ent researchers [3|-|6| . Low coercivity is desired in the core 
materials (or it is only used in recording for a short life- 
time) , which means low hysteresis loss per cycle of opera- 
tion. Contrary to this, high coercivity is used in applica- 
tions where magnetic recording media is frequently used 
or need to have a long life. This constant strategy is the 
main trend-setter in magnetic thin film design especially 
in recording applications |7|. In other respects, control- 
ling thickness of film to obtain the magnetic hysteresis 
at a right shape to suit desired technological purpose has 
been the key issue in disguise 0. Therefore, the core 
meaning of how the spontaneous spin formation mecha- 
nisms involved in these reduced structures are different 
form the bulk have become the focus of frequent inves- 
tigation issues. However, to the best of our knowledge, 
there is not any systematic investigation regarding the 
hysteresis behavior and its infiuences on each layer of the 
ferromagnetic thin film due to the underlying complexity 
of the reduced dimension. So, there are some important 
issues which remain in suspense; frequency dispersion of 
HLA with corresponding remanence and coercivity for 
different layer indices in parallel with the classification of 
dispersion curves, topological evolution of the hysteretic 
response of each layer with respect to the modified sur- 
face exchange, etc. 

Consequently, the aim of the study is to give a more 
complete picture of the overall hysteretic behavior of fer- 
romagnetic thin films especially at different values of 
modified exchange interaction for a fixed temperature 
value (the temperature has been selected a certain cri- 
terion). For this purpose, we first investigate how the 
hysteretic properties of a film, including their relevant 
dispersion curves, depend on the dynamic system param- 
eters by means of effective-field theory (EFT) . After that 
we construct various kinds of hysteresis (in fact, they are 
collections of curves well-selected from the frequency dis- 



persion of the HLA) to confirm the results and to observe 
how the response depends on the layer index. Finally we 
give our results by discussing the magnetization reversal 
mechanism within reasonable bounds. 



II. METHODOLOGY 

From the theoretical point of view, one of the mod- 
els more widely used to study the magnetic properties of 
surfaces is semi-infinite Ising model due to the crossover 
of dimensionality and their strong uniaxial anisotropy, 
so the system can be modeled by a layered Ising-type 
structure which consists of interacting L parallel layers. 
Each layer is defined as a regular lattice with coordina- 
tion number z — A. The familiar Hamiltonian of the thin 
film is given by, 
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where Si is the spin operator at a lattice site i and any 
spin variable can take the values Sj = ±1. As is known, 
(. . .) subscript bracket symbolizes the nearest neighbor- 
ing in first summation. The second summation is over all 
the lattice sites. The exchange interaction Jij between 
the spins on the sites i and j takes the values according to 
the positions of the nearest neighbor spins. Two surfaces 
of the film have the intralayer coupling Ji. The inter- 
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layer coupling between the surface and its adjacent layer 
(i.e. layers 1, 2 and L, L — 1) is denoted by J2- For the 
rest of the layers, the inter layer and the intralayer cou- 
plings are assumed as J3 . The system has three exchange 
interactions where Ji, J2, Ji > favors a ferromagnetic 
alignment of the adjacent sites as shown in Fig. (1) which 
has been presented in our prior work (see Ref. [81]) and 
the Zeeman term describes interaction of the spins with 
the field of the sinusoidal form 



h{t) = /iQ cos(a;i). 



(2) 



where t is the time and h^ is the amplitude of the oscil- 
latory magnetic field with an angular frequency w. 

Our system is in contact with an isothermal heat bath 
at given temperature T. So, the dynamical evolution 
of the system may be given by non-equilibrium Glauber 
dynamics [9J based on a master equation. Number of L 
different representative magnetization time series for the 
system can be given by usual dynamical EFT equations 
which are obtained by differential operator technique |10|, 
lll| . In order to handle the multi-spin correlations, a 
decoupling approximation (DA) Iz] can be used as 
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..sf^) = {sf^)...{s'.P)...{sf^)^ (3) 



which is essentially identical to the Zernike approxima- 
tion 13] in the bulk problem, and it has been successfully 
applied to a great number of magnetic systems including 
the surface problems [l4lla |. Thus, the dynamical equa- 
tions of motion for each layer are in the form of 
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Here ruk which is the magnetization time series of fcth 
layer is defined as 



TOfc = (S 

with the coefficients 
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1, 



,,L, 



An = cosh(J„V)/(a;)|:r=o, 
En = sinh(J„V)/(x)|:r=o, 



(5) 
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where n — 1,2,3 and V = d/dx is one dimensional dif- 
ferential operator and the function f{x) is given by 



f{x) = tanh[/3(a; + h{t))] 



(7) 



Here j3 — l/ksT and ks represents the Boltzmann con- 
stant. The effect of the differential operator on an arbi- 
trary function /(x). 



exp(aV)/(a;)]:c=o = f{x + a)\^^o, 



(8) 



with any real constant a. perbolic trigonometric functions in terms of the exponen- 

tial functions we get the most compact form of Eq. ^ 
By using the binomial expansion and writing the hy- as 
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Here, l/r is the transition per unit time in a Glauber 
type stochastic process and throughout our calculations 
it is selected as t = 1 for simplicity. Self-consistent non- 
linear ordinary differential equations in Eq. ([9]) are solved 
by fourth order Runge-Kutta method (RK4) to get the 
evolution of the mk{t) by regarding the equation as an 
initial value problem. 

The system has three dependent Hamiltonian vari- 
ables, namely frequency of external magnetic field cj, am- 
plitude ho and the film thickness L. For certain values 
of these parameters, temperature and the Ji, J2, J3 in- 
teraction constants, RK4 will give convergency behav- 
ior after some iterations i.e. the solutions have property 
m{t) = m{t + 2tt/uj) for arbitrary initial value for the 
magnetization {m(t = 0)). Each iteration, i.e. the calcu- 
lation of magnetization at i -f 1 from previous magneti- 
zation for t, is now performed for these purpose whereby 
the RK4 iterative equation is being utilized to determine 
the magnetization for every i. In order to keep the itera- 
tion procedure stable in our simulations, we have chosen 
10^ points for each RK4 step. Thus, after obtaining the 
convergent region and some transient steps (which de- 



pends on Hamiltonian parameters and the temperature) 
the layers' average magnetization can be calculated from 



}k^ -^ j> mk{t)dt 



(11) 



where rrik is a stable and periodic function and finally 
the dynamical order parameter (DOP) can be calculated 
by the arithmetic mean as 



1 ^ 

Q = 7 X Q* 



(12) 
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As cut-off condition for numerical self-consistency, we de- 
fined a tolerance 



Qk\t-2Tv/u ~ Qk\t 



t + 2lT/L0 



<10' 



(13) 



meaning that the maximum error as difference between 
the each consecutive iteration should be lower than 10~^ 
for all step. On the other hand the HLA corresponding 
to energy loss due to the hysteresis is defined as 



mk{t)dh{t) — hoU} f mk{t)sm{ujt), (14) 



^(fc) 



and the modified surface exchange interaction has been 
defined to determine the different characteristic behavior 
of the system in certain range as 



Ji = J3(l-I-A,). 



(15) 



There are three possible states for the system, namely 
F, P and the coexistence phase (F+P). The total magne- 
tization time-series ■m{t) in convergent region is satisfied 
by this condition 



m{t) = ~m{t + tt/uj) 



(16) 



in the P phase which is called the symmetric solution. 
On the other hand, in the F phase, the solution does 
not satisfy Eq. (fTH]) and this solution is called as non- 
symmetric solution which oscillates around a non-zero 
magnetization value, and does not follow the external 
magnetic field i.e. the value of Q is different from zero. 
In these two cases, the observed behavior of magnetiza- 
tion is regardless of the choice of initial value of magne- 
tization m(0) whereas the last phase has magnetization 
solutions symmetric or non-symmetric depending on the 
choice of the initial value of magnetization corresponding 
to the coexistence region where F and P phases overlap. 
The main goal of the treatment is adjunctly manifesting 
the frequency dispersion of the HLA and corresponding 
coercivity (He ) with remanence {Mr ) of layers. Note 
that, more details about formulation and the reasons for 
selection procedure of the related values of system param- 
eters at which the calculations carried out can be seen in 
Ref. i. 



III. RESULTS AND DISCUSSION 

With the guidance of the work presented by Akta§ et 
al. Q, the frequency dispersion of the HLA and the 
corresponding coercivity with remanent magnetization 
for a pure crystalline ferromagnetic thin film with thick- 
ness L = 7 are globally presented in Figs. dH), ^ and 
([3]) . As a result of evaluations by the conviction that we 
have reached, the dependency characteristics with vary- 
ing frequency has an intersection at a critical value of 
the modified exchange interaction. For instance. Fig. ([T]) 
shows the hierarchical reversal of the HLA in each layer as 
changing A^. At a fixed ho, it can be seen that the surface 
has the lowest loop area in A^ < A* till the symmetry 
loss occurs as well as the innermost layer has the lowest 
loop area in the opposite regime (namely A^ > A*). In 
other respects, the layer indices are of no importance in 
frequency dependency at the critical value of modified 
exchange interaction, i.e. if the modified exchange inter- 
action ever had it's critical value all the layers would have 
the same hysteretic behavior. Areal magnitude sequence 
of hysteresis loop for the layers becomes reversed between 
each other at A*. Physical mechanism mentioned above 
can be briefiy explained as follows: At a fixed /iq and 
As value in As < A* regime but on increasing layer 
index k (e.g. k — 1 corresponds to the surface layer) 
hysteresis loop becomes larger at low w but it slightly 
gets smaller at high w. This is since for inner layers 
there are relatively more neighboring per magnetic-sites 
which cause locally larger magnetic interaction. With a 
stronger magnetic interaction in that local lattice region. 



it becomes more difficult for spins to follow the field and 
the hysteresis almost tends to be asymmetric. Surface 
atoms are embedded in an environment of lower symme- 
try than that of the inner atoms and consequently the 
exchange constants between atoms in the surface region 
may differ from the bulk one for experimental evidence. 
The lower As regime corresponds to surface type of mag- 
netic ordering with this aspect. However, one can reverse 
the whole aforementioned scenario by controlling (actu- 
ally by increasing) As, i.e. in As > A* regime increas- 
ing layer index k hysteresis loop becomes smaller at low 
cj but slightly larger at high uj. This is since for inner 
layers there are relatively more neighboring but far less 
exchange constants per magnetic-sites which cause rela- 
tively smaller magnetic interaction than that of the sur- 
face one anymore. Finally, theoretical evidence forced us 
to remark that the free surface couldn't break the trans- 
lational symmetry since the hysteretic and/or magnetic 
properties of the free surfaces seriously overlap with the 
bulk one at the critical value of As (namely A*). The 
investigations depicted in Figs. dH) and ^ are also con- 
sistent with the whole arguments about the mechanism. 

At a fixed ho and As, for purely symmetric region, with 
increasing w, Ak increases at low lu whereas it decreases 
as w whereas it decreases in agreement with the well- 
known behavior in literature where frequency at the peak 
ujc corresponds to the phase-lag. In As < A* regime, this 
Wc shifts to lower values of cu as increasing the layer index 
k due to the locally stronger magnetic interaction. With 
the same reason, at a fixed uj, A^ increases (especially at 
low oj) with increasing layer index k for As < A* while 
reverse of the overall arguments given above are valid for 
As > A* regime. One can easily observe the intersection 
on ujc values of layers via changing As parameter. 

In addition, in Fig. ^ and both Fig. © with (O, the 
coercivity He increases with increasing the field ampli- 
tude ho whereas increasing in remanence Mr can be 
seen in significantly asymmetric region. It means that 
HLA Ak increases with increasing ho and this increment 
is dependent of the other system parameters. This is 
since the higher amplitude provides higher magnetic en- 
ergy to the system and this higher energy provides more 
magnetic force in causing the spins to follow the field. 
Consequently, the frequency at the peak (wc) in symmet- 
ric region for a selected layer index and value of As shifts 
to higher uj as the phase-lag between ordinary magnetiza- 
tion and field signal is smaller for the higher amplitude. 

Next, the hysteresis loop collection of each layer for re- 
lated system parameters and selected the best three ap- 
propriate values of frequency are presented in Fig. ([¥]). 
As is acquainted with the low frequency limit, the mag- 
netization can be more saturated and this shows itself 
as the dominant magnetization processes in the different 
regions of the curve (i.e. reversible boundary displace- 
ments, irreversible boundary displacements and magne- 
tization rotation). The figure shows that the hysteretic 
response of each layer clearly has the same character- 
istics at the critical value of the modified exchange A* 
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FIG. 1. (Color online) The frequency dispersion of HLA for 
each layer of the film with thickness L = 7 at a fixed tem- 
perature fcsT/Js = 3.5. The curves as seen in each group 
of triple-column have been plotted for the selected amplitude 
representatives as /lo = 1-0 and 2.5. They profiles across 
the film differ qualitatively in two regimes (As < A* and 
As > As). Each inset is stated in terms of visual clarity. 



independently with the other parameters. Hierarchical 
reversal of the hysteresis loop of each layer for various 
values of amplitudes and/or frequency also can be seen 
in Fig. (Ul by following the group of three columns in- 
cluding the related regimes of Ag ( A^ < A* , Ag = A* 
and As > A*). One can observe the effects of symme- 
try on this reversal mechanism (following ability of the 
layers to external field sweep) by comparing the repre- 
sentative amplitudes between each other. At the fixed 
As = A* value, for h^ — 2.5 representative, it can be 
seen that the hysteresis loop has a saturated s-shape and 
tends to increase in size with increasing w at low fre- 
quency region e.g. u < 0.3. However, on further increas- 
ing w, the loop gets its maximum area and later reduces 
to an oval shape with its major axis parallel to the field 
axis. This is the result of phase-lag between the ordinary 
magnetization and the field signal. At very low oj, the 



FIG. 2. (Color online) The frequency dispersion of coercive 
field for the same selected system parameters of the film with 
thickness L — 7 depicted in Fig. ((T}. Each inset is stated in 
terms of visual clarity. 



field period is large and the magnetic spins have sufficient 
time to follow the field signal so the phase-lag between 
the ordinary magnetization and field signal is small. The 
system is strong enough to follow the external perturba- 
tion at those limitations, hence the hysteresis loop looks 
like a slim s-shape. However, on increasing w, the field 
sweeps faster and the spins have less time in following the 
field. Consequently, the phase-lag gets larger and so does 
the hysteresis. At a particular point, when the phase- 
lag approaching its critical value, the hysteresis gets its 
maximum size. After that, if lu still increases in small 
doses, the spins feel very difficult in following the field so 
the overall magnetization does not change much as the 
field changes. Therefore the hysteresis turns its shape 
in oval-shape with its major axis lying along the field 
axis. Additionally, for Hq = 1.0 representative, a rare 
pre-symmetry loss morphology can be seen by following 
the topological evolution of hysteresis loop with increas- 
ing frequency at the fixed Ag > A*. For instance, as seen 
at id = 0.07 value, an exotic left-handed s-shape beyond 
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FIG. 3. (Color online) The frequency dispersion of remanent 
magnetization for the same selected system parameters of the 
film with thickness L = 7 depicted in Fig. ([1]). 



the oval form temporally takes place between maximal 
area frequency Wc and the frequency value of symmetry 
loss which is absolutely stated about uj > tUc- The result 
of general assessment about the hysteresis collection as 
follows: For low frequencies, the hysteresis loop is in the 
form of saturated s-shape and tends to increase in size 
with increasing uj at low frequency region. At a critical 
LUc value, the loop gets its maximum area. On further 
increasing lo, hysteresis enters into the process of smooth 
directional- veering and it has a left-handed s-shape until 
the symmetry loss. Forthwith, the symmetry loss occurs 
at a characteristic frequency value and after that point, 
the hysteresis has an asymmetric shape. 

The process stated between at the beginning of 
directional-veering and symmetry loss deserves a 
scrutiny. As seen in Figs. ([T]), ^ and (jS]), symmetry 
loss takes place at a gibbous-like eminentia. This phe- 
nomenon has been zoomed in Figs. ([T|) and ([2]) for better 
visibility. Particularly, as seen in Fig. ([3]), both of the 
positive and negative remanence appears at low frequen- 
cies, but only positive values survive at high frequencies 



due to the symmetry loss. The related eminentia shows 
itself more specifically on right arm of the mutual co- 
ercivity and it is far more distinct in surfaces for low 
amplitudes in condition with A^ > A*. Notice that the 
encountered symmetry loss characteristics for all obser- 
vations on this problem are the member of Type ILb 
family of frequency dispersion presented by Akta§ et al. 
(iTJ I . Low amplitude relatively provides smaller magnetic 
energy to the system and this weak energy could not 
provide enough magnetic force in causing the spins to 
follow the field. This means the surfaces fall all over it- 
self to adopt to the asymmetric phase in the A^ > A* 
regime, so, the surfaces have extremely left-handed sat- 
urated s-shape hysteretic response within that period. 
It is also of interest to consider the topological evolu- 
tion of the hysteresis as changing uj to investigate how it 
relates with spin formation (namely magnetic domains) 
growth and/or nucleation mechanism. Two branches of 
the hysteresis curve intersect each other at loc and the 
HLA has the global maxima related to the multiplying 
of its two-axis. Therefore, the magnetic moment remains 
approximately constant value in the half-time period of 
external field. At the fixed values oi uj < uJc and uj > uJc 
respectively, corresponding hysteresis curves are topolog- 
ically contra-oriented mutually and completely different 
growth mechanism is onset. According to times-series 
of the magnetization versus the external magnetic field, 
increasing the field frequency at first, obstructs the satu- 
ration of the ordinary magnetization due to the decreas- 
ing energy coming from the oscillating magnetic field in 
a half-time period which facilitates the late stage do- 
main growth by tending to align the moments in its di- 
rection (i.e. the magnetization begins to fail following 
the oscillatory field) and this makes the occurrence of 
the frequency increasing route to continuous dynamical 
phase transition due to the incomplete reversal of mag- 
netic moments. The aforementioned directional-veering 
mechanism in the A^ < A* regime begins to lose its 
significance. Since as the A^ is lowered more and more 
then the dipole-dipole interaction-induced energy contri- 
bution gets smaller. The quantity of the energy contri- 
bution which comes from the dipole-dipole interaction 
corresponds to the self-volition to remain in the dynami- 
cally symmetric phase of the system. Hence, the surfaces 
can stay in the symmetric phase in the smaller frequency 
until the frequency of the external field becomes greater 
than the relaxation time of the system. 
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FIG. 4. (Color online) The related HL collection of each layer for the film with thickness L — 7. Modified exchange interaction 
As ranges from 0.0053 (first column) to 0.3053 (second column) and 0.6053 (third column) with varying frequency as uj = 0.07, 
0.3 and 0.53 for related amplitude representatives at ksT/Js = 3.5. 



IV. CONCLUSION 



As a complementary work, the effect of As on the topo- 
logical evolution and strength of HLA for each layer has 
been discussed by means of effective field theory based 
on a standard decoupling approximation for dynamical 
fsing-type thin films of L = 7 driven by an external os- 
cillatory magnetic field with inner coordination number 
z = 4. The time evolution of the system has been pre- 
sented by utilizing a Glauber type stochastic process. For 
pre-selected values of the system parameters, our investi- 



gation has been focused on frequency dispersion of HLA, 
coercivity and the remanence of each layer. 

Concretely speaking, in order to propound a link be- 
tween domain growth and/or nucleation mechanism and 
symmetry loss in magnetic thin films, we appeal to hys- 
teresis for its topological evolution in respect to the fre- 
quency for throughout our work. Hence, some unpreten- 
tious improvements on the relation has been developed 
in accordance with the explanations given before in lit- 
erature by some experts. 

In conclusion, we hope that the results obtained in this 
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work would shed light on the further investigations of the 
dynamic nature of the critical phenomena in pure crys- 
talline ferromagnetic thin films and would be beneficial 
from both theoretical and experimental points of view. 
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